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I. INTRODUCTION 
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We investigate the evolution of the density-density correlations in the isoscalar critical condensate 
formed at the QCD critical point. The initial equilibrium state of the system is characterized by a 
fractal measure determining the distribution of isoscalar particles (sigmas) in configuration space. 
Non-equilibrium dynamics is induced through a sudden symmetry breaking leading gradually to 
the deformation of the initial fractal geometry. After constructing an ensemble of configurations 
describing the initial state of the isoscalar field we solve the equations of motion and show that 
remnants of the critical state and the associated fractal geometry survive for time scales larger than 
the time needed for the mass of the isoscalar particles to reach the two-pion threshold. This result 
is more transparent in an event-by-event analysis of the phenomenon. Thus, we conclude that the 
initial fractal properties can eventually be transferred to the observable pion-sector through the 
f^*) ■ decay of the sigmas even in the case of a quench. 
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Experiments of a new generation, with relativistic nuclei, at RHIC and SPS are currently under consideration with 
^vq . the aim to intensify the search for the existence and location of the QCD critical point in the phase diagram of 
;> strongly interacting matter Important developments in lattice QCD Q and studies of hadronic matter at high 

temperatures [4] suggest that the QCD critical endpoint is likely to be located within reach at SPS energies. It is 
therefore desirable to explore the range of baryon number chemical potential [Ib =100-500 MeV by studying collisions 
at relatively low energies, s/s ~ 5-20 GeV per nucleon pair. A decisive observation, in these experiments, associated 
with the development of a second-order phase transition at the critical endpoint of QCD matter is the establishment 
of power laws in momentum space (self similarity) in close analogy to the phenomenon of critical opalescence in QED 
matter [5J. These power laws reflect the fractal geometry of real space, and a characteristic index of the critical 
behavior is the fractal mass dimension Df which measures the strength of the order-parameter fluctuations within 
the universality class of critical QCD [6[ . 

The physics underlying the endpoint singularity in the QCD phase diagram is associated with the phenomenon of 
Qh. chiral phase transition, a fundamental property of strong interactions in the limit of zero quark masses. In this case 
and for a given chemical potential [ib there exists a critical temperature T cr above which chiral symmetry is restored 
and as temperature decreases below T cr , the system enters into the chirally broken phase of observable hadrons. It is 
believed that, for two flavors and zero quark masses, there is a first-order phase transition line on the (T, (1b) surface 
at large \ib |Z|- This line ends at a tricritical point beyond which the phase transitions become of second order. In 
the case of real QCD with non-zero quark masses, chiral symmetry is broken explicitly and the first-order line ends 
at a critical point beyond which the second-order transitions are replaced by analytical crossovers [§[ . 

The order parameter of the chiral phase transition is the chiral field $ = [a, 7?) formed by the scalar, isoscalar 
field a together with the pseudoscalar, isovector field tt = (ir + , ir°, 7r~). Both fields are massless at the tricritical 
point and when the symmetry is restored at high temperatures, their expectation value vanishes, (7?) = (cr) = 0. 
However, in the presence of an explicit symmetry breaking mechanism (non-zero quark masses) the pion and sigma 
fields are disentangled at the level of the order parameter of the QCD critical point, which is now formed by the sigma 
field alone. The expectation value of the cr-field remains small but not zero near the critical point so that the chiral 
symmetry is never completely restored. The valuable observables in this case are associated with the fluctuations of 
the cr-field, {5a) 2 ~ (c 2 ), and they incorporate, in principle, the singular behavior of baryon-number susceptibility 
and in particular the power-law behavior of the cr-field correlator Q . 

The aim of this work is to study the evolution of critical correlations during the development of the collision, 
in the neighborhood of the QCD critical point. In the initial state we assume that the system has reached the 
critical point in thermal equilibrium so the cr-field fluctuations are described by the 3-D Ising critical exponent 
5 and in particular by the fractal dimension Df = (5 ~ 5). The crucial question from the observational point 
of view is whether in the freeze-out regime, which follows the equilibration stage, the relaxation time-scale r re i of 
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these fluctuations is long enough compared to the time-scale Tth associated with the development of a massive cr-field 
beyond the two-pion threshold (m a > 2m w ). Both time-scales (T re ;,r t ^) are characteristic parameters of the out-of- 
equilibrium phenomena (cr rescattering) which take place during the evolution of the system (towards freeze-out) and 
the requirement r re i 3> Tth guarantees that critical fluctuations may become observable in the cr-mode (cr — > tt + tt~ 

0-IHJ)- 

In order to quantify these effects, we adopt in this work the picture of a rapid expansion (quench) which is a realistic 
possibility in the framework of heavy-ion collisions. We study the out-of-equilibrium evolution of the initial fractal 
characteristics of the cr-field and we search for time scales r re ; , r t h satisfying the above constraint in a particular class of 
events (event- by- event search). The dynamics of the system is fixed by a two-field Lagrangian, £(cr, ff), together with 
appropriate initial conditions. The out-of-equilibrium phenomena are generated by the exchange of energy between 
the cr-field and the environment which consists of massive pions initially in thermal equilibrium. 

In section II the formulation of the problem and in particular the equations of motion, the initial conditions and 
the generation of thermal 7r-configurations are presented. In section III numerical solutions of the evolution of critical 
fluctuations are given and discussed whereas in section IV our final results and conclusions are summarized. 



II. FORMULATION OF a FIELD DYNAMICS 



In our approach we assume an initial critical state of the system in thermal equilibrium, disturbed by a two- 
field potential V(a, tt), in an effective description inspired by the chiral theory of strong interactions [111 hj. The 
3-dimensional Lagrangian density is 

C = hid^a + d^d^) ~ V(a, 7?) (1) 

with the potential 

Til ^ 

V(a, vr) = — (a 2 + tt 2 - v 2 ) 2 + ^(n 2 - 2v Q a + 2v 2 ) , (2) 

where a — cr(x, t) and 7? = tt(x, t). The potential has the usual tr-model form plus a term which breaks the symmetry 
along the tr-direction. With the addition of the mass term for the pion field, we ensure that it has a constant mass 
equal to m n . Finally, the constant terms in shift to zero the value of the potential at the minimum. We fix 
the parameters of the Lagrangian using the phenomenological values m^w 139 MeV and vo ~ 87.4 MeV, w hereas 
the known uncertainty in the phenomenological value of zero temperature cr-mass, given by m a = ^ l /2X 2 v^, yields 
10 < A 2 < 20 for 400 < m a < 600 MeV. 

The equations of motion resulting from ([I]) arc: 

<7 - V 2 fT + A 2 (er 2 + tt 2 - Vq)<t = v ml 
#-V 2 7f + AV 2 +7? 2 -^)7? + m 2 7? = 0, (3) 

where tt 2 = (tt+) 2 + (tt ) 2 + (tt") 2 . 

Using a constant value vq in eq. [2] implies a non-vanishing mass (finite correlation length) for the cr-field, m a — 
y // 2X 2 VQ, in contradiction with its critical profile. This inconsistency is restored if we assume a finite-time mechanism 
instead of the instant quench i.e the instantaneous formation of the potential ([2]) at t — 0. The simplest model which 
leaves the equations of motion unaffected, used also in cosmological phase transitions [H|, is the so-called linear 
quench [I3 |. It assumes that the minimum v of the potential increases linearly with time, starting from zero and 
ending at the zero temperature value vq w 87.4 MeV after a time interval r which is the quench duration: 

v(t)=v t/r for t<r 



v(t) = w 



for 



t > T. 



(4) 



In this way we acquire m CT = at t — 0, as expected for the critical cr-field and, with increasing time, m a approaches 
its zero temperature value. 

To proceed numerically we have to discretize eqs. j3j on a lattice. We use the following leap-frog discretization 
scheme: 
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and similarly for the other three equations considering the 7r-field. In eq. ([5]) a is the lattice spacing while dt is the time 
step. The upper indices indicate the time instants and the lower indices the lattice sites. As usual we perform an initial 
fourth order Runge-Kutta step to make our algorithm self-starting, and we impose periodic boundary conditions. 

We are interested in studying the evolution of the above system using initial field configurations dictated by the 
onset of the critical behavior. In this case we expect that the cr-field, being the order parameter, will possess critical 
fluctuations, and the 7r-fields to be thermal, while the entire system will be in thermodynamical and chemical equilib- 
rium. Obviously, the subsequent evolution, determined by eqs. ([3]), will generate strong deviations from equilibrium. 
Before going on with the detailed study of the dynamics, we first describe in the following subsections the generation of 
an ensemble of field configurations on a 3-D lattice possessing the characteristics of the critical system. This ensemble 
enters in the subsequent analysis as the initial condition. 



A. Generation of initial ensemble of critical a-configurations 

The absolute value of the cr-field introduced in the previous subsection is interpreted as local density, and the 
corresponding critical behavior is described by a fractal measure demonstrated in the dependence of the mean " mass" 
M(xq, R) on the distance R around a point Xq defined by: 



M(x ,R) = { (6) 

JR 

obeying the power law 

M(xq, R) ~ R Df (7) 
for every xq. Df is the fractal mass dimension of the system (HEEEli 

and the mean value is taken with respect to 
the ensemble of the initial cr-configurations. The production of the a configurations building up the critical ensemble, 
characterized by the fractal measure given in eqs. (|6I7|) , has been accomplished in [l8|. In fact the fractal properties 
of the critical system can be produced as an ensemble average, through the partition function: 

Z = j 5[a]e-^\ (8) 

with r[er] the scale invariant effective action at T = T cr , fi = fi cr : 

T[a} = Jj° x {l(Vaf + ga s+1 }. (9) 

The sum over field configurations in (jHJ) can be saturated through saddle point solutions of T[cr], which consist 
approximately of piecewise constant configurations extended over domains of variable size. A detailed description of 
the corresponding simulation algorithm is given in (l9| . Using this algorithm, after averaging, the mean "mass" in 
the ensemble satisfies eq. with a fractal mass dimension 

For the 3-D Ising universality class, D = 3, the isothermal critical exponent is S « 5, and the coupling g ss 2 [20| . 
therefore Df « 5/2. 

The power-law behavior of M(Sq, R) = (J R \a(x)a(xo)\ d 3 x) around a random xo, averaged inside clusters of volume 
V, is illustrated in fig.[T] The M(xq, R) versus R figure is drawn as follows: For a given xo of a specific configuration 
we find R of the cluster in which it belongs, taken ~ \rV, and we calculate the integral / \a(x)a(xo)\ d 3 x, thus 
acquiring one point in the M(xo, R) vs R figure. For the same xo we repeat this procedure until we cover the whole 
ensemble, and the aforementioned figure is formed. Averaging in xq leads only to a slight modification of the result, 
since M(xq, R) ~ M(xq + l,R), with I spanning the entire lattice, therefore in the following we replace M(xo, R) by 
M(R). We observe that in the log- log plot of M(R) vs R, the slope 4", i.e the fractal mass dimension Df according 
to ([7]), deviates from the theoretically expected value of 5/2 by less than 1%. Furthermore, the mean value of the 
a- field (spatial average) is almost zero as expected to happen near the critical point. 

With this procedure we acquire an ensemble of field configurations and the expected power law arises as a statistical 
property after ensemble averaging [I9j ]. Alternatively, one could extend the notion of the fractal dimension using 
individual configurations, where the quantity of interest is the exponent df of the power law of the integral 

d 

\U {X)U \XQ ) \ U X U X'O ~ JrL 

R. 



m(R)= / \a(x)a(x )\d A xd J x - R af . (11) 
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FIG. 1: M(R) — (J R \a(x)a(xo)\ d 3 x) versus R for the ensemble of a -field configurations. The slope 9 , i.e the fractal mass 
dimension Dj is equal to 5/2 within an error of less than 1%. 



Calculating df for each configuration we obtain a distribution around 5/2 with standard deviation ss 0.05. As expected 
the ensemble average of df is (df) = Df, within an error of less than 0.5%. 



B. Generation of initial thermal ^-configurations 

We generalize the method of (2lj | in order to produce an ensemble of 3-D 7r-configurations in real space, corresponding 
to an ideal gas at temperature Tq. The unperturbed Hamiltonian for the classical scalar field theory in three dimensions 
is 

H=\J d 3 x[(d t 7r(x, t)) 2 + (Vtt(x, t)) 2 + ml-K{x, t) 2 ]. (12) 
The free particle solutions for t — are 



7r(a;,0) = / — 7Tfc e = / — -== — e 



(2tt)3 ^ u J_ x (2tt)3 ^%T k 

^■<» = / + °°7^&oe^= [ +X ^M i{a *_ k -a k )e^. (13) 



(2tt)3 ™ J_ x (2tt)3V 2 



where uj k = \/k 2 + rri*. 

Now, choosing an initial classical density distribution [2l[ 

p[*,w] =Z- 1 (J3 ) exp {-#, H[n, *■]}, 

and substituting the Hamiltonian (fT2")) with the free particle solutions (TT3]), we finally get 

f / ,+ °° d 3 fc 1 
p[x k ,y k ] = Z~ 1 {fi ) exp|-/3 y -—^(^ + yjg) j, ( 14 ) 

with /3o = 1/Tq, and where: a k = x k + iy k with x k ,y k real. In order to produce a thermal ensemble (at temperature 
To) of configurations for n(x, 0) and 7r (af, ) , we select and from the gaussian distribution (fT4")l . assemble afc 
and then substitute in (|13p . Lastly, since we have three components of the pion pseudoscalar field, we independently 
repeat this procedure accordingly. All the characteristics of the 7r-ensemble such as the correlation function (7r(x)7r(x + 
5x)) — (tt(x)) (n(x + 53)}, which turns out to be a 5-function, are consistent with the assumption of an ideal thermal 
gas. 
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III. NUMERICAL SOLUTIONS 

We study the evolution of the system determined by equations ([3]) which we solve in 3-D 25 x 25 x 25 lattice, using 
as initial conditions an ensemble of 10 4 independent cr-conhgurations on the lattice generated as described above, 
i.e possessing fractal characteristics, and 10 4 configurations for each it component corresponding to an ideal gas at 
temperature To « 140 MeV. The initial time derivatives of the cr-field, forming the kinetic energy, are assumed to be 
zero, since this is a strong requirement of the initial equilibrium. The used population is by far satisfactory since the 
results converge for ensembles with more than 7 x 10 3 configurations (numerically tested) for the considered lattices 
(sizes from 15 x 15 x 15 to 25 x 25 x 25 sites). In addition, we find that the obtained results are independent of 
the lattice spacing a despite the discontinuities in the derivatives of the piecewise constant configurations, as the 
corresponding variation 5a goes to zero fast enough so that the limit lim Q ^o ^ exists [l9[ . Lastly, varying the value 
of To between 100 MeV and 180 MeV it has a negligible effect on the results. We investigate the evolution of M(R) of 
the whole ensemble which initially follows a power law ~ with \I/(0) = Df = 5/2. We also study the evolution 

of m(R) for each configuration. At t — m(R) possess a power-law behavior of the form ~ R^(°\ with the exponent 
ip(0) = df normally distributed around 5/2 (with standard deviation ss 0.05). 

We evolve the system for various A 2 , corresponding to different m a values at T = 0, and for various quench times 
r. In figs. [2] and [3] we demonstrate the evolution of the mean field values (a) and (ir + ) (the other components are 
similar), where the averages are taken over all statistically independent configurations. As expected (a) oscillates 
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FIG. 2: Time evolution of the mean field value (a) for various X 2 and r cases. 



around the potential minimum which relaxes after r, while (tt + ) stays around zero due to the absence of a linear in 
7? term in the potential J2]). 

In fig. [4] we present M(R) vs R for three successive times and we observe that the slope ^(t) fluctuates, leading to 
a distortion of the initial fractal geometry. In fig. [5] we depict the evolution of the slope ^f(t) of M(R) versus R (each 
^f(t) value obtained trough a linear fit), for the same A 2 and r values as before. With the dashed line we plot the 
time average of ^(t) defined by = j J ^(t')dt' . We observe the remarkable phenomenon that the characteristic 
exponent after reaching the value of the embedding dimension 3, it fluctuates and for particular times it becomes 
almost equal to \& (0) = 5/2. Thus, after the first deformation, the initial critical behavior of the whole ensemble is 
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FIG. 3: Time evolution of the mean field value (n) for various X 2 and r cases. 



t=0fm 



t=0.25 fm 



t=0.5 fm 



>P(0)=D f =2.52 



>F=2.92 



¥=2.78 



R(fm) 



FIG. 4: M(R) vs R for three successive times, t = fm, t — 0.25 fm and t — 0.5 fm, for typical X 2 and r values. The solid 
lines mark the power-law fit, and the dotted line has slope 2.5 and is used as a reference. 



partially restored and deformed repeatedly. A detailed explanation of this revival is given in [22j |. The key point is 
that the partial restoration takes place when (<r(t)) passes through its lower turning point, where the cr-field (seen as 
a system of coupled oscillators) reaches a state similar to the initial one. This phenomenon weakens gradually, and 
finally the dynamics dilutes completely the initial critical behavior. 

This phenomenon is also visible in the evolution of the slope ip(t) of m(R) versus R, for each configuration. Indeed, 
in fig. [6] we demonstrate the evolution of the mean field value as well as of ip and of its time average, for A 2 = 20 
and r = 1 fm, for three independent configurations of the ensemble, corresponding to initial value ^(0) = df equal 
to 2.51, 2.50 and 2.53 respectively. We observe an oscillatory behavior of ip(t) similar to that of 'f(i) of the whole 
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FIG. 5: Time evolution of (solid line) and (^)t (dashed line) for the whole ensemble, for the same A 2 
before. The vertical line marks the time t = Tth where m a reaches the threshold 2m^. 
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ensemble presented in fig. [3 

In a heavy ion collision, if the fireball passes near the critical point we expect a to have almost zero mass. As the 
system expands, its temperature decreases, v increases towards its zero temperature value Vq and m a = y / 2\ 2 v(t) 2 
tends to its freeze out value. However, when it reaches the threshold m a — 2m T , at time t — r t h, it starts decaying 

into pions, and the decay rate of a — > ir + ir~ , being proportional to the available phase space parameter \ — 

becomes larger as m a increases. So if the critical characteristics of a have survived at that threshold, they can be 
transferred to the produced pions, leaving signatures of the critical point at the detectors. Note here that these pions 
are not affected by the thermal ones which constitute the environment. In figs. [5] and [6] the vertical line depicts the 
threshold time t = Tth- We are interested in investigating the time averaged measures after this threshold since these 
quantities are observable. 

The first measure one has to look at is the evolution of ^(t) and (^>}t for the whole ensemble, after t = r t h, as it is 
shown in fig. As we observe these vary between 2.6 and 2.9, for the various A 2 and r values, offering weak traces of 
the initial power law of 5/2. This cogitation is amplified by the fact that if we evolve our system with conventional 
initial conditions (Df « 3), then \T/(t) and (\l/)t remain always w 3, as we have tested, supporting the assumptions 
that slopes < 2.9 indicate initial critical behavior. However, slopes close to 3 could originate from conventional strong 
processes leading to power- law correlations with very small intermittency exponents [23| . Therefore we have to refer 
to more sophisticated measures in order to acquire a direct observation of the QCD isothermal critical exponent 5, 
independently of the specific model. One way to achieve this is to perform event-by-event analysis of our results. 

As we have mentioned, we use an ensemble of a- field configurations, each one possessing its own df which can 
vary, leading to a distribution with mean value « 5/2 and standard deviation w 0.05. An event-by-event analysis 
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FIG. 6: (a), ip(i) (solid line) and (tjj)t (dashed line) for three independent configurations, for A 2 
vertical line marks the time t = Tth inhere m a reaches the threshold 2m 7T . 
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of the system evolution consists in calculating the percentage P of the initial configurations (all of which have 
tp(0) = df rj 5/2) that possess again this value at time t (actually we count those with 2.4 < tp(t) < 2.6). Obviously, 
P(t) leads to experimentally accessible effects after t = r t h- In fig. [7] we depict the evolution of P(t), and its time 
average (P)t for t > Tth (the vertical line marks t = Tth), for the same A 2 and r values described above. As we observe, 
P(t) is quite large and decreases with time as expected, while it presents some fluctuations which correspond to the 
peaks of ^(t) evolution of fig. [5] As we mentioned above we can define the relaxation time-scale T re i, after which the 
initial critical behavior is completely lost, as the time where (P)t becomes less than 0.5%. Therefore, if r re i is larger 
than Tth then the critical characteristics will be transferred to the produced, through the decay a — > t + tt~ , pions. 

The quench duration seems to affect the results significantly. For quench time r = 1 fm the aforementioned time 
averaged (P)t can be quite large for the first 1-2 fm after t = Tth, as can be observed in fig. [JJ and it reaches the cut 
0.5% at T re i « 40 fm, a clearly larger value than the expected freeze-out time in heavy ion collisions. In other words, 
T re i ~ 40 fm means that (P)t > 0.5% at the freezeout. For r = 2 fm, (P) t has a similar behavior but r re i decreases 
to « 1.6 fm, for the A 2 = 20 case. For r = 3 fm, (P) t is always below 0.5% since in this case t = r t h is reached 
after the first large fluctuation of ^(t), as it is shown in the e and f plots of fig. \5\ Thus, in this case, according to 
the above definition r re i = 0. In general for faster quenches, (P)t, i.c the time average of the percentage of the field 
configurations that reacquire the value ip(t) w 5/2 after t = Tth, is significantly larger. On the other hand, for slower 
quenches Tth increases and it supplants r re ;, thus the dynamics dilute the traces of the initial critical profile, before 
the decay a — > 7r + 7r _ activates. The dependence of the ratio T re i/Tth on the quench duration r is depicted in fig. [5] 
for the A 2 = 20 case. As we observe, for a wide range of r values (r < 2.1 fm in the A 2 = 20 case), r re ;/r t ^ > 1 and 
especially for small r it can be sufficiently large (T re i/T t h — 10 for r ~ 1.6 fm). As the model parameter r increases, 
the ratio T re i / Tth decreases almost exponentially as can be induced from the exponential fit presented with the dashed 
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FIG. 7: Time evolution of the percentage P(t) of the configurations that reacquire the value ip(t) ~ if>(0) = df « 5/2 oi t. With 
the dashed line we depict its time average (P) t after t = Tth- 



line in fig. EJ 

The effect of A 2 (that is the value of the m a at T = 0) is not so crucial and T re i shows a small decrease with 
decreasing A 2 , while Tth increases slightly. However, for a given value of the interval r re i — Tth the pion production 
from the critical ex's is controlled by A 2 . Since the c-decay rate is larger for larger ma/m^, we expect faster decay 
for A 2 = 20 compared to A 2 = 10, since in the former case the final m a is larger. Experimentally we expect a larger 
number of produced pions in the A 2 = 20 case. Inversely, when A 2 = 20 a smaller interval r re ; — r t h is needed in order 
to produce a given number of pions through the decay of critical er's. This is the main difference of the two A 2 cases, 
although the corresponding plots look alike. 

The discussion above reveals that the critical behavior of the cx-field will sustain, in a significant percentage of the 
initial configurations, for times after the cr's start to decay, thus transferring the critical profile to the produced pions, 
despite the abrupt non-equilibrium evolution from its initial equilibrated state. 

It is important to exclude the inverse possibility, i.e appearance of events with fractal geometry being generated out 
of the evolution of conventional initial conditions. Therefore, we perform the following test: We evolve our system 
using an ensemble of ^-configurations with random initial conditions. This conventional profile always possesses 
(df) w 3 with a very small standard deviation (w 0.02). In this case not even one configuration (out of 10 4 ) acquires 
slope < 2.96 ever, while the distribution becomes even narrower as times passes. As a consequence, P(t) and (P)t are 
always exactly equal to zero, within the used population. Therefore, the detection of events with critical profile at 
the freeze out, will offer safe signatures of the initial criticality of the system. The requirement (P)t > 0.5% remains 
robust for a long time and for a variety of A 2 and r as we have shown, providing a guarantee that such events are 
very likely to be observed in experiments of relatively high statistics. 

The above treatment shows, together with other approaches in the past @, 0, HH HH , that a systematic study 
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FIG. 8: The dependence of the ratio T re i/rth on the quench duration r for A 2 = 20. The dashed line depicts an exponential fit. 



of event-by-event fluctuations in relativistic nuclear collisions gives rise to a powerful tool in the search for novel, 
unconventional effects associated with QCD phase transitions (critical fluctuations, disoriented chiral condensates, 
dissipation at the chiral phase transition). In this work we have examined the evolution of critical fluctuations 
neglecting the effects of dissipation in this out-of-equilibrium process. Although a detailed study of dissipative 
aspects of this evolution is beyond the scope of this work, an estimate of their influence on the fractal structure of 
critical fluctuations is necessary before drawing our final conclusions. For this purpose we have considered a simplified 
modification of a field dynamics adding a dissipative term Ta(x, t) — jv(t) in the first equation ([3]). v(t) is the noise 
term fulfilling the relations (v(t)) = and (v{t)v{t')) — j 2 S(t — t'). The coefficients T and 7 are related through the 
formula 7 = y/2TT/V, resulting from the fluctuation-dissipation theorem, with T the environment temperature and 
V the system volume. The friction coefficient T introduces a new time scale in the process (rt = associated 
with the rolling down of the effective potential [25[ . With this modification we have solved the equations of motion 
Q following the method discussed in the beginning of this section. 

We use the typical values A 2 = 20 and r = 1 fm, and for the friction we impose T" 1 ~ 0.5 fm, which corresponds 
to the strong friction regime [2^]. As we observe in fig.EJ an d ^(t) change significantly compared to the non- 

dissipative case, and the traces of the initial fractal geometry dilute faster. However, the change in (P) t does not spoil 
its observability since Tth and r re i are almost unaffected (see the inset in fig. Eh). Hence, dissipation influences the 
mean values evolution, but it does not affect qualitatively the fluctuation-related measures such as (P)t- Therefore, 
the signals of the initial critical characteristics are still visible in the event-by-event analysis. Moreover, we have found, 
as expected, that weakening the dissipation, i.e decreasing T, its effects are suppressed and the system evolution tends 
smoothly to that of the non-dissipative case. 



IV. SUMMARY AND CONCLUSIONS 



In this work we have explored the dynamics of critical fluctuations which are expected to occur in the c-mode, near 
the QCD critical point. For this purpose we have adapted the a- model Lagrangian in order to describe correctly the 
characteristics of the order parameter associated with the critical endpoint of the QCD phase transition [12] ■ The 
issue is of primary importance in the search for the existence and location of the QCD critical point, in experiments 
with nuclei. At the phenomenological level these fluctuations are expressed through the fractal mass dimension of the 
(T-field configurations, determining the properties of the condensate at criticality [6(. We have studied the evolution 
of the initial critical characteristics of the tr-field in thermal equilibrium and the possibility to reveal signals of critical 
fluctuations not affected by the dynamics which drives the system out-of-equilibrium. We found that for a wide range 
of quench time-scales (r) and coupling values (A 2 ), the initial critical profiles may survive for a long time after the 
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FIG. 9: (a), ip(t) (solid line), (ip)t (dotted line) and P(t) (solid line), (P)t (dotted line), in the strong dissipative case (V 1 ~ 0.5 
fm), for X 2 — 20 and r = 1 fm. 



system reaches the two-pion threshold value of the sigma mass, m a — 2m v . This result is more transparent in an 
event-by-event study of the phenomenon, where the evolution of individual configurations is investigated. 

As a consequence of this study, the fractal dimension of critical fluctuations in QCD matter, turns out to be a 
remarkable index for the location of the QCD critical point, in experiments with nuclei, of relatively high statistics. 
In fact, it remains robust, in a class of events, against dynamical effects and it leads to a characteristic pattern of 
intermittent fluctuations 0, [23[ in transverse momentum space, providing us with a signature of the critical point 
without ambiguities owing to dynamics. 
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